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The proton structure function F2 is studied in the low x regime using BFKL evolution. The next to 
leading logarithmic (NLL) analysis requires the inclusion of running coupling effects which lead to 
off-diagonal terms in the evolution kernel. An all-orders resummation is used to improve the collinear 
behavior of the NLL BFKL result. We emphasize the theoretical uncertainties that appear throughout 
the analysis and give a comparison to the combined HERA data. 

1 Introduction 

In 2010 HERA made public the combined results (TJ obtained by HI (2) and ZEUS (3J Collaborations 
for the proton structure function F 2 (x, Q 2 ) at low values of the Bjorken x variable and a quite broad 
range of values of the photon virtuality Q 2 . This observable became specially convenient to test the 
region of applicability of the theory based on the high energy or Regge limit, which corresponds to 
the center of mass energy of the system y/s being asymptotically larger than any other scale involved. 
In Deep Inelastic Scattering (DIS) a hard scale is provided by the high virtuality of the photon. Since 
the x variable is given within very good approximation by the ratio between the photon virtuality and 
the center of mass energy squared, we can refer to the Regge limit as the limit of low Bjorken x. In 
this regime large logarithms of energy appear, dominating the scattering amplitude, and they need to be 
resummed to all orders. Such a resummation is achieved by the so-called BFKL evolution equation |4]- 

The aim of the present study is to analyze the theoretical uncertainties encountered in the determina- 
tion of F 2 at NLL accuracy using high energy factorization [9 |. Care has to be taken when introducing 
running coupling effects and it is also needed to resum to all orders the leading collinear singularities 
which are numerically large in this kinematical region flO[jT4[ . Concerning the running of the coupling, 
we compare the results obtained using a running with a Landau pole with a model which freezes in the 
infrared and is compatible with power corrections to jet observables [15]. A numerical analysis of the 
gluon Green function in the kinematic region of interest is provided. We conclude with a preliminary 
comparison of our theoretical calculation for F 2 with HERA data. 
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2 Analyzing the proton structure function 
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High energy factorization allows to write the proton structure function as 
a convolution in transverse momentum space of a nonperturbative object 
describing the proton (proton impact factor $ p ) with the photon (photon 
impact factor $ 7 ), calculated using perturbation theory, together with a 
gluon Green function /, linking both process-dependent components and 
accounting for the BFKL evolution: 
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Fig. [T] shows the different parts involved. Although an analytic expression 
for the photon impact factor at next to leading order accuracy (16 ■ 18 1 has 
been recently calculated fl9) we use for simplicity in our analysis the lead- 
Figme 1 m g or( j er resu it as presented in p0| . The proton impact factor needs to be 

modeled. Our results are based on a simple choice which introduces just a few parameters to be deter- 
mined by the experiment and it presents a Poissonian-like distribution in transverse momentum space 
with its maximum around the confinement scale. Finally, the gluon Green function is governed by the 
BFKL equation. Its LL solution is smooth and convergent but not sufficient to explain the DIS data. 
The first attempt to have a good description of F2 would consist on studying the next order of accuracy. 
However it is known that the NLL kernel is unstable in collinear regions. We have found that by intro- 
ducing an all-orders collinear resummation consistent with the NLL solution following the procedure 
given in 1 14 1 we not only eliminate the collinear instabilities but also get a good preliminary description 
of the data. Figs.[2]and[3]compare the LL gluon Green function to the collinearly improved one. Since 
this formalism does not modify the NLL results but only gives higher order corrections there is certain 
freedom in the way of performing the resummation. In the present analysis we use an expression for the 
NLL eigenvalues which includes the action on the impact factors of the differential operators breaking 
the scale invariance of the kernel. 

At NLL accuracy one needs to account for the running 
of the coupling. This gives an analytical expression for the 
kernel which contains a differential operator in the Mellin 
variable 7 pT) . There is in principle no theoretical re- 
striction (other than having an hermitian hamiltonian) on 
whether to act with this operator in a symmetric way [22 1, 
i.e. on both proton and photon impact factors or in an asym- 
metric way (only on one of them). Nevertheless, it turns out 
that each option produce very different results, as shown in 
fig. [2] The reason for this is that each of them naturally leads 
to a different scale for the running coupling due to the scale 
of the logarithm accompanying the term in (3q of the kernel, 
responsible for the running. This is, however, an assump- 
tion about higher order terms again. We could have decided 
to leave the logarithms without absorbing them into the ex- 




Figure 2: Action of the differential op- 
erator. Upper curve: LL; mid curve: 
asymmetric differential operator acting 
on the photon; and lower curve: symmet- 



ric choice. 

pression for the running coupling. We have also compared the results obtained for two different models 
of the running, a perturbative one, with a Landau pole, and the one presented in 1 15 1 and described 
earlier in the introduction. However, as it can be seen in fig. 3b we are studying a region in which these 
are minor effects, since we ask the transverse scale to be perturbative. 
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(a) Dependence on k 2 for sq = Q 2 (solid lines) and 
so = kq (dashed lines) at LL (upper set, green) and NLL 
collinear improved (lower set, red) accuracy. 
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(b) Dependence on the model for the running: the solid 
lines correspond to the IR finite running and the dotted ones 
to the perturbative one at LL (upper set, green) and NLL 
collinear improved (lower set, red) accuracy. 



Figure 3: Numerical analysis of the gluon Green function. 



The last important point of discussion in this analysis is what to do with the choice of energy scale 
so appearing in the gluon Green function: 
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It is known that any dependence of the cross section on this scale must cancel at NLL accuracy. However, 
if we want to express / as a function of a; a shift in cj is produced leading to a remaining dependence on 
it that appear as higher order corrections. A natural choice in this case would be the DIS scale, so = Q 2 , 
so that (s/sq) u = The symmetric choice as a product of the internal scales s — kq was used to 



calculate the NLL BFKL solution [21]. As it can be seen in fig. 3a there is a difference in the results 
obtained with each version. 

Figure [4] shows one of the possible preliminary fits that can be done of Fi within this theoretical 
framework. In particular, a symmetric version of the differential operator has been used together with 
the symmetric energy scale so = kq, the perturbative running coupling and proton impact factor given 
by 

* p (k 2 )^(k 2 /Q^ e kVQ^ 

The expression for the photon impact factor used can be found in J23| . 
A detailed analysis of the work here presented can be found in [24|. 



Acknowledgements 

The European Comission (LHCPhenoNet PITN-GA-20 10-264564) is acknowledged for funding the 
expenses of the conference. 



DIS 2012 



3 



us 

1.4 
1.2 
2 1.0 
0.8 
0.6 

1 x 10" 4 2x 10" 4 5x 10" 4 0-001 0.002 0.005 0.010 

Figure 4: Preliminary fit to F 2 with 6 = 1.246, Q 2 = 0.368GeV 2 and A p = 0.07346. 
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